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ABSTRACT 


An accurate method has been developed for predicting effects 
of streamline curvature and coordinate, system rotation on tur- 
bulent boundary layers. A new two-equation model of turbulence 
has been developed' which serves' as the basis of the study. In 
developing the new model, physical reasoning is combined with 
singular perturbation methods to develop a rational, physically- 
based set of equations which are, on the one hand, as accurate 
as mixing-length theory for equilibrium boundary layers and, 
on the other hand, suitable for computing effects of curvature 
and rotation. The' equations are solved numerically for several 
boundary layer flows over plane and curved surfaces. For incom- 
pressible boundary layers, results of the computations are gen- 
erally within 10^ of corresponding experimental data. Somewhat 
larger discrepancies are noted for compressible applications, 
although the overall level of accuracy is obscured by uncertain- 
ties in the experimental data. Results of the study lend further 
confidence to the notion that curved streamline effects can be 
accurately predicted with second-order closure of the' turbulent- 
flow equations of motion. 
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DEFINITION 

Constant In the law of the wall 
Skin friction 

Parameters in turbulence model equations 
Specific turbulent mixing energy 
Nondimenslonal mixing energy 
Nondlmenslonal streamfunction (Eq.l20) 
Roughness height 
Reference length 
Mach number 
Static pressure 

Production term in dissipation-rate 
equation 

Radius of curvature 

Reynolds number based on plate length, 
momentum thickness, displacement thickness 
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Static temperature 

Vdloclty component in x direction 

Boundary-layer-edge velocity, freestream 
velocity 

Velocity component in y direction 
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DEFINITION 

Nondlmensional dissipation rate 
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Coordinate normal to a solid body 

Parameters in turbulence model equations 
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solution 
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Momentum thickness 
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Parameter in turbulence model equations 
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Fluid mass density 
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Similarity parameter in defect-layer 
solution 


vi 

f 


V 



NOTATION (camelu'de'd) 


SYMBOL 

4 > 


to 


DEFINITION 
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solution 


Specific turbulent dissipation rate; 
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1. INTRODUCTION 


In a prior study, Wilcox and Chambers^ found that the magnitude 
of streamline curvature effects on turbulent boundary layers 
can be predicted reasonably well with second order closure of 
the turbulent flow equations of motion. The study showed that 
streamline curvature and coordinate-system rotation effects can 
be understood on a fairly fundamental basis only by considering 
the interactions of the various components of the Reynolds stress 
tensor and the mean flow. Hence, most of the effort focused upon 
developing a turbulence model based on the Reynolds stress equa- 
tion (RSE). The most Important result of the Wilcox-Chambers 
study is the straightforward physical Interpretation of stream- 
line-curvature effects inferred from the analysis. As a brief 
review, the RSE model predicts that for boundary-layer flows the 
following things are true: 

1. In equilibrium the Reynolds shear stress is 
proportional to the normal Reynolds stress perpen- 
dicular to the plane of shear, <v’^>. 

2. <v*^> is strongly affected by streamline curvature 
and system rotation. 

3. The large changes in <v*^> occur mainly because of 
RSE transformation properties and pressure-strain 
correlation terms. 

While other researchers have suggested that <v’^> plays a key 
role in curved-surface boundary layers, no one previously demon- 
strated that curvature effects can be accurately predicted by 
accounting for changes in <v*^> caused by curvature. 
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As a corollary result of the study, understanding of curvature/ 

rotation effects gained from analyzing the RSE model led to 

2 — j . 

Improvement of the simpler Saffman ttw.oeqq.uat.i'ontfeu'hhai'e.nee 

model. Because of Its Inherent simplicity (relative to the RSE 

model), the two-equation model appears to have the best potential 

for serving as an engineering design tool. 

While these developments are most encouraging, results of numer- 
ical applications of the RSE model and the Saffman two-equation 
model leave something to be desired. For example. Figure 1 cora- 
pare§ computed and measured' velocity profiles for two equil- 
ibrium boundary layers with adverse pressure gradients. . Both 
the RSE and Saffman models yield Inflected profiles In the 
defect-layer portion of the boundary layer. Since for equil- 
ibrium boundary layers the RSE model reduces almost Identically 
to the Saffman model. It Is unsurprising that both models yield 
the same curious Inflections. With the models performing so 
poorly for equilibrium boundary layers. It Is difficult to gen- 
erate much confidence In model predictions for such complicated 
effects as those of curvature and rotation. In the present 
study, application of the RSE model to flow over a segmented 
spinning body (see Appendix) further demonstrates the RSE model's 
deficiencies In the defect layer. 

The objective of the present study has been to Identify and 
eliminate the source of Inaccuracy In the Saffman model. Since 
the Saffman model serves as the foundation upon which the RSE 
model is built, elimination of deficiencies in the Saffman model 
should lead to Improvement of the RSE model. The study consists 
of two segments. The first segment (Section 2), focuses primar- 
ily upon two-equatlon-model defect-layer structure and equil- 
ibrium (plane-surface) boundary-layer applications; a new baseline 
two-equation model is developed. Then, In Section 3 , with suit- 
able curvature-rotation modifications the baseline model equations 
are solved numerically for several boundary-layer flows over curved 
surfaces. Section 4 discusses study results. 
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Figure 1. Comparison of computed and measured velocity profiles for turbulent 
boundary layers with adverse pressure gradient. 





2. DEVELOPING A BASELINE MODEL 


This section traces the evolution of a new two-equation model 
of turbulence. First, Saf fman-model-predlcted defect-layer 
structure for a flat-plate boundary layer (FPBL) is analyzed 
using slngular-perturbat ion ■ techniques . Next, a new model is 
postulated and is shown to yield defect-layer properties which 
are in much closer agreement with measurements than those of 
the Saffman model. The high Reynolds-number form of the 
Saffman model, the new model, and two other popular two-equation 
models are then compared by computing several boundary layers. 

f 

Finally, with appropriate viscous modifications, the new model 
is applied to equilibrium boundary-layer flows. 

2.1 PERTURBATION ANALYSIS OF THE DEFECT LAYER 

As discussed in the Introduction, the problems in Saffman and 
RSE model solutions apparently originate in the defect layer. 

In order to determine the nature of the models* deficiencies, 
a singular-perturbation analysis of Saffman-model-predicted 
defect-layer structure has been conducted for a zero pres- 
sure gradient FPBL. As a first step in the analysis, the 
. Saffman equations are suitably scaled and simplified for 
the defect layer. Then numerical solutions are presented in 
which the boundary-layer edge conditions are parametrically 
varied. Finally, it is found that the Saffman model predicts 
a too-large peak value of the turbulent length scale, £. 

2.1.1 Formal Expansion Procedure 

The Saffman model equations of motion for incompressible two- 
dimensional boundary layers are 

Viscous modifications are changes required to add molecular 
viscosity to the momentum and model equations, where appropriate 
to permit integrating the equations through the viscous sublayer 
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where u and v are velocity components in the strearawlse, x, and 
surface-normal j y, directions , respectively; is the boundary- 
layer-edge velocity; v is kinematic viscosity and the eddy 
viscosity is the ratio of turbulent mixing energy, e, and 
turbulent dissipation rate, w. The six empirical parameters 
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a, a*, g, g*, a, a* are given by 

a = a* = 0,50 
g = .0,15 , g* = 
a* =0.30 

a = g/a*-4aK^ = O.I638 for k = 0,.4l 
where < is Karman’s constant. 



Equations 1-5 must be solved subject to boundary conditions at 
the boundary-layer edge (y^°°) and at the surface (y-»-0). In the 
defect layer, we anticipate that the eddy viscosity is propor- 
tional to Ug<S* where 6* is displacement thickness, wherefore 

V ^ V _ 1 /fr 

e/w Ug6* Reg* . . 
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where Re^* denotes Reynolds number based on 6*. Since Reg* >> 1 
for turbulent boundary layers, molecular viscosity will be neg- 
ligible in the defect layer. Consequently, defect-layer solu- 
tions must be subjected to boundary conditions as y-»-0 which are 

valid at the outer edge of the viscous sublayer, namely, the law 

.fe ' 

of the wall. As shown by Saffman and Wilcox, this means 
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Equations 7 are valid on the sublayer scale as u^y/v^~. 
Formally we are matching the outer (defect-layer) solution 
to the inner (sublayer) solution. 

Finally, at the boundary-layer edge we impose the values of 
u, e and w, viz: 


u = 
e = 
w = 



as y " 


( 8 ) 


Following Bush and Fendell?^ (wet&nte§^ueey|l3ecfpilgwing';non^ 
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where L is a length typical of the distance from the plate 
leading edge. Note that the eddy viscosity becomes 


e 
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U 6* ^ 
e W 


( 10 ) 


and that the quantity U^6*/u^ is the Clauser thickness. It is 
convenient to work with the streamf unction, ijj, defined by 


iL = 11 


U 3 3y 
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Finally, we Introduce the scaled streamfunction, F, given by 

^ ^ U~ ( 12 ) 

e 

In terms of the transformation defined in Equations 9-12, the 
equations of motion become (neglecting molecular Viscosity): 
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where the quantity m Is defined by 
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To solve these equations, we assume that the velocity deviates 
only slightly from the freestream velocity which Is expressed 


by writing a perturbation expansion In powers of u /U << 1 as 

X ^ 

follows : 
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Similarly, the scaling for E and W follows from Equation 17 so 
that we write 
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Then Introducing 


Ui = 9F^/3n 


( 21 ) 


and assuming u^/U^ ~ 6*/L, the momentum equation becomes, to 
leading order 
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Equation 22 will have self-similar solutions (i.e., independent 
of 5) provided the following four quantities are Independent 
of ? : 
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where c„ is skin friction, x is shear stress at the wall, and 
f w 

p is pressure. Substituting Equations 23-26 into Equation 22 
yields 
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As shown by Bush and Pendell, 

p 

Wy = p(l) as Rsg* ■»■ “ 

while the momentum integral equation shows that if self-similar 
solutions exist, and 3,p are related as. follows: 


ct,j, = .1 + 3B,p 


Hence, assuming that ahd 3,p are constant, the defect-layer 
equations are: 
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2.1.2 Numerical Solution of the Defect-Layer Equations 

Equations 30-34 were solved using an Implicit, second-order 
accurate, tlme-marchlng, finite-difference code. Only the zero- 
pressure-gradient case (Bip=0) was considered. To achieve con- 
verged solutions, we found It necessary to Impose the boundary- 
layer-edge conditions at a sharp turbulent-nonturbulent Inter- 
face.^ This Is unsurprising as Saffman ' has shown that the model 
equations predict the existence of such Interfaces when e and m 
approach zero In the freestream. The Interface structure remains 
unaltered when e and w have finite values In the freestream. That 


+ 

The sharpness manifests Itself as a discontinuity In 3u/9y and 
. 3u)/9y. 
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is, as a generalization of the analysis by Saffman, upon approach 
ing a turbulent-nonturbulent interface (y=6) from the turbulent 
side, the solution behavior is 

e-e^ (6-y)^^ 

asy->6 

u-Ue (6-y) ) 

provided a = a* = 1/2. 

Several computations were done to determine solution sensitivity 
to e^, 0 ) , a, and a*. Figure 2. shows velocity profiles for a 
set of computations in which solution sensitivity to e^ was 
examined. In the computations, the following experimental 
boundary condition, was Implemented in an attempt to eliminate 
solution dependence upon oj^ : 


|| = 0 . at y = 6 (35) 

where 

Z = e^^^/u ■ (36) 

As shown in the figure, the model predicts only a very weak 
"wake" component. The computed velocity profile deviates sig- 
nificantly from Coles' composite "wall-wake" profile. Values 
of a*e /u^ less than 0101 (corresponding to a turbulence Inten- 
slty of about 0.3^) have little effect on the solution. Larger 

values of a*e /u^ cause the solution to deviate from Coles' 
ex 

profile even more than the e^^O solutions. Since Coles' pro- 
file provides a fairly accurate representation of experimental 
data, the model apparently is quite Inaccurate in the defect 
layer. 
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Figure 2. Comparison of computed Saffman-model velocity profiles with 
Coles* composite wall-wake profile. 







No combination of a, a*, e^, and oj^ was found which would give 
more accurate defect-layer velocity profiles. The source of 
the difficulty can be seen by Inspection of the computed length- 
scale, ly profile. In terms of the eddy viscosity, s = e/w, 
can be written as 

E = e^/^il (37) 

As shown by Ng and Spalding^ (who use Equation 31) » J<// a*' 
should be nearly equal to the standard mixing length provided 
e is close to the turbulent energy, as it is for the Saffman- 
model solution. Figure % compares the computed length-scale^ 
profile with the mixing-length profile deduced by Klebanof f . 

The peak value of the mixing length is between 0.095 and 0.106. 
By contrast, the peak value, of £// predicted by the Saffman 
model is about 0.1767 


Further numerical experimentation showed that the peak value 
of H/6 can be suppressed somewhat by using extremely large 
values of This is generally accompanied by severe numer- 

ical difficulties and by little Improvement in the velocity 
profiles. A second means for reducing the peak value of Jl/6 
is to reduce the value of a; the reduction is slight, however, 

The value of a* has virtually no effect on (Z/S) although 

irio.x 

the value of 6 depends strongly upon a*. 


Hence, the Saffman model inherently predicts too large of a 

value of (A/6) in the defect layer. Consequently, since e 
ma.x 

is computed with reasonable accuracy, the eddy viscosity is 
overestimated as can be seen by noting that the Saffman model 
typically predicts 


^max 

U 6* 
e 


'x- 0.035 


(38). 
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a value more than double. that assumed (0.0168) in the mlxing- 
length/eddy-viscosity formulation. The large peak values of 
1/6 and e/U^6* predicted by the Saffman model cause the curious 
inflected velocity profiles we have so often obtained with the 
model (see Figure 1). Consistent with the defect-layer analysis 
presented above, accurate skin friction, shape factor, etc., can 
only be obtained by using unreasonably large values of co at the 
boundary-layer edge with an attendant distortion of the velocity 
profiles . 


2.2 A NEW TWO-EQUATION TURBULENCE MODEL 

After obtaining numerical solutions to the defect layer equations 
and discovering the exaggerated peak value of i-/6 , the next logi- 
cal step was to find the specific cause of the Inaccuracy. In 

0 7 

order to do this, we examined the Jones-Launder*^ and Ng-Spalding 
turbulence models, both of which predict more-realistic peak values 
of a/6. As shown by Wilcox and Chambers,^ these two models can be 

written in terms of the Saffman variables (e,w) as follows: 
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tion term by P , 
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that r ■l^<< (if) ignore the term proportional 
6 oy dy 

Jones-Launder model, and (§) denote the produc- 
we can write the following single equation which 
three (including Saffman's) models: 



(41) 


Saffman model 
Jones-Launder model 
Ng-Spaldlng model (42) 

es-Launder and Ng-Spalding 
models have an additional dissipation term proportional to 
0 )^O£,/ 3 y)^. This term has the effect of augmenting dissipa- 
tion near a surface where 3Jl/9y is largest; by contrast, in the 
dofect-layer where 3il/3y is smallest, the net dissipation of 
will be smaller than near the surface. The effect of this term 
is hence consistent with the notion that the most intense dis- 
sipation occurs in the smallest eddies which appear nearest a 
solid boundary. In the defect layer only large eddies are 
present wherefore dissipation will be smaller than near the 
surface. 


where 
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Equations 4l and 42 show that the Jon 


Because the Jones-Launder and Ng-Spalding models predict more- 
realistic peak values of 2/6, the possibility exists that the 
Saffman-model might be Improved if such a term is added to the 
equation. A value for x can be chosen by noting that, in 
the wall layer 

(y) " ^^3) 
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Hence with c^= 702, k= 0.41, a* = 0.3, a=l, we obtain 


6 - = 6 
^ (3il/8y.)^ 


c a*2 
w 

a 


2.43 


(44.) 


It thus appears that x should lie between about 2.0 and 2.5. 
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0 = 0 *= 1/2 
3 = 0.15 , e* = 0.09 
. a* = 0.3 

a = 3/a* - (4-x)cr.K^ 


(47) 


Figure 4 presents the computed Jl/6 profile for x=2* The peak 
value of S./6 Is much closer to the Kiebanoff mixing-length data 
than that obtained with the Saffman model. Computed velocity 
profiles are presented in Figure 5- As shown, excellent agree- 
ment between the numerical profiles and Coles* wall-wake profile 
is obtained. Additionally, the mixing energy is in somewhat 
closer agreement with the measured RMS vertical velocity fluc- 
tuation data of Klebanoff^*^ (see Figure 6). Using X“2*5 causes 
only slight changes in computed profiles. 
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Figure 4. Comparison of computed length-scale 
profiles with mixing-length profiles 
deduced from measurements. 
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Figure 5. Comparison of computed Basellne-modelr velocity profiles 
with Coles' composite wall-wake . prof lie . 





Figure Comparison of computed turbulent 

mixing energy profiles with Klebanoff's 
RMS vertical velocity fluctuation data. 



It is possible to obtain skin friction from the perturbation . 
solution. Matching demands that as n-^0, the defect-layer solu- 
tion must equal the limiting form of the sublayer solution as 
u^y/v-^oo. Now, Inspection of the velocity profiles in Figures 2 
and 3 . shows that 
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- ^ log n + a] 


as 


(il8) 


where A. is a constant which can only be determined as part of 
the solution (note that Bush and Pendell mistakenly demand that 
A=0). In terms of the sublayer variables. Equation 48 can be 
rewritten as 


U^ K ^ V 
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, (49) 


Finally, from the sublayer solution. 
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where B-5.5 is the standard constant in the law of the wall. 
Matching yields the following condition: 


^ = B + A + ^-logReg* (51) 

T 

wherefore the relation between skin friction, c^ = 2u^/U| and 
Reg* becomes 


2k‘ 


(log Reg*)2 + 
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(A+B.)^ 
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Figure 7 compares skin friction for 3-^^ X"2 with the Hopkins- 
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Inouye correlation, in which we have assumed the shape factor 
to be 1.30 wherefore 
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Figure "S- Comparison of computed skin friction with the Hopklns-Inouye 
correlation of experimental data. 




^ = 3.221 (log Re^*)^-H- 9.215 log Reg* + 3.373; 

Hopkins-Inouye (53) 

As shown In the figure, .the Saff man-model is generally about 

10/S higher than the Hopkins-Inouye value. This is consistent 
with previously computed values of c^. The B'asi$irierainddel’>^'wlth 
X=2 predicts values of c^ much closer to the Hopkins-Inouye 
values. It hence appears that the proposed revision holds 
promise for significantly improving Saf fman-model accuracy. 

2.3 TWO-EQUATION MODEL COMPARISON 

To further test the Baseline model, a series of boundary-layer 
computations were performed. Since the objective of this study 
was to establish the most solid foundation possible for future 
research with the . RSE model, we left the possibility open that 
a formulation other than the Baseline or Saffraan models might 
be more accurate. Hence, a comparison of the Baseline, Saffman, 
Jones-Launder , and Ng-Spalding turbulence models was conducted. 
There is no need to consider additional models as all current 
two-equation models are variants of these four. 

This comparison provides an objective basis for determining 

which of these models is most accurate overall. While previous 

12 13 

comparisons of two-equation formulations have been made, ’ 
no one has made more than a cursory examination. Until now, no 
one has presented parallel numerical solutions for the various 
models. In this study, the models were tested under exactly the 
same conditions using the same numerics, boundary conditions, 
and starting profiles . Hence from the results of these compu- 
tations, it is possible to determine precisely how well the 
models perform relative to each other and which model should 
serve as the basis for further analysis of streamline curvature/ 
system rotation effects on turbulent boundary layers. Table 1 
lists the flows considered and the data sources used. 
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Table 1. Plows Computed in Two-Equation- 
Model Comparisons 


PLOW 

DATA SOURCES 

• Plat Plate Boundary Layer 

Hopklns-Inouye^^ skin fricti.dn 


correlation; Klebanof f flat 


plate dataj Wieghardt^^ flat 


plate data. , 


1 ^ -.14 

• Bradshaw Adverse Pressure 

Bradshaw data; Coles’ ver- 

• Gradient 

slon of Bradshaw data. 

• Aridersen Adverse Pressure 

4' 

Andersen' data. 

Gradient 


• Ludwieg-Tillmann Favorable 

14 

Ludwieg-Tillmann data; 

Pressure Gradient 

Coles 'l4 version of Ludwieg- 


Tillmann data. 

• So-Mellor Constant Pressure 

So-Mellor^^ data; Meroney^^ 

Plow over convex wall 

correlation . 


Before proceeding to results of the computations. Subsections 2.3-1 
through 2.3.3 present (a) the equations which constitute the vari- 
ous models, (b) boundary conditions used in the computations, and 
(c) Initial profiles. 

2 . 3.1 The Model Equations 

The inviscld-incompresslble forms (l.e., negligible Mach riumber 
and molecular viscosity) of the model equations were solved. The 
inviscid (l.e., high Reynolds number) forms were chosen because 
of, the nonuniversality in the way viscous effects are Included 
in the various turbulence, models . In this way we help eliminate 
any preferential treatment of a given model. For all models, 
the conservation of mass and momentum equations are identical, 
viz. 
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Mass Conservation 


9u _j_ ^ 
3x 3y 


0 


( 54 ) 


Momentum Conservation 


9u , 

"K * 



1 dp ^ r 3u ~ 
p dx 3y 3y. 


( 55 ) 


The two model equations for each formulation are as follows. 


Saffman Model 


Turbulent Mixing Energy 


3e , 3e 




= e*“|e + 


( 56 ) 


Dissipation Rate 


8a) . 3o)‘ 

+ ''W 


= " ihtl ‘57) 


The values of the constants are 

3 =.15 3 * =.09 a = 3 /a*- 4 cK^ = . 1638 - 

a = . 50 a* = . 50 a* = . 30 

while the ratio of e and o) is the kinematic eddy viscosity, e, 
i. e . , 


e = e/o) 


( 58 ) 


The Baseline model differs from the Saffman model only by the 
addition of one term to the dissipation rate equation as 
explained in the previous section. 
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Baseline Model 


Turbulent Mixing Energy • 


( 59 ) 


Dissipation Rate 

9W^, 9u)^ r |9U| d 1 2i 9 r 9(0 ^*] f ^ ^ \ 3 / Cr\\ 

‘ [“1^1 ■ H" ^r.3y J ■ ^'’(sy) “ 


Note that x been set equal to 2. All constants are the same 
as in the Saffman model except for a which is now given by the. 
following expression: 

a = B/a* - 2 o.k^ 


The kinematic eddy viscosity is given by Equation 58. 


As noted earlier, similar to the Saffman and Baseline models, the 
Jones-Launder model uses turbulent mixing energy for one "turbu- 
lence density"; the second turbulence density is the dissipation 
function, e^. 


Jones-Launder 


Turbulent Mixing Energy 


,.9e , ,. 9e _ _/^9u \ 
^9x ^9y V3yj 


£ + -i-/— 

d 9yVag3yj 


(61) 


Dissipation Function 


u- 


^_fd 

9x 


9e 


+ V- 


d _ ^<fd [m\^ §d^. . 9 /e ^Sd ^ 

9y ^ i / ' "'"52ee”'’ ^y \g^ ;9y 


( 62 ) 
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where 


c^=1.55, C 2 = 2, . a^=l,0, a^=1.3 

and. the kinematic eddy viscosity takes on the following form: 

e = .09 (63) 


Finally, the Ng-Spalding model uses turbulent mixing energy and 
turbulent length scale. 

Ng-Spaldlng 

Turbulent Mixing Ene’i^'gy 


9e , 


ie ^ c ■ (64) 

9y . \3y/ D ^ ^ 


Length Scale 


u~(eJO + v^(e£) = Cpj2,e 


2 

/9u\ _ 

\9y/ 


c + c (-V )e^^^ + 
i m w\y/ ) 


_3_r^ 

3yL^2 


ay 


(eS,) 


where 


Cj^= .09, c^= .98, c„ = .059, c,,, = 702, 0^=1. d, a^=1.0 


m 


w 


and 


£ = e ^'2 ^ 


(66) 


All four models use a turbulent mixing energy equation. Com- 
monly, e is assumed to be the turbulent kinetic energy. 

However, Wilcox and Chambers have shown^ that more appropriately 
e should be interpreted as 9/4<v’^> where v' is the fluctuating 
component of the velocity component normal to the surface. The 


(65) 
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above fact is essential when determining such things as 
starting profiles. 

The second turbulence parameter differs for the Saffman and 
Baseline (m), Ng-Spalding (A), and Jones-Launder (e^j) models. 
The eddy viscosity, e, is used in the same manner in each model 
as an integral part of approximating the Reynolds shear stress. 
Hence the various expressions for e can be equated and rela- 
tions between w, and Jt can be derived. It is necessary to 

have such relations in order to impose the same boundary layer 
edge boundary conditions and initial profiles for each of the 
models. In terms of e and the dissipation rate and dissipa- 
tion function are 

CO (0= =e (67 ) 

and 

e^f:= (68) 

In addition to comparing the four models, results of the compu- 
tations yield two important conclusions about all the models. 
One, accurate starting profiles can be obtained by knowing only 
the mean velocity profile and the Reynolds shear stress profile 
at the starting location. Two, the models are not sensitive to 
either initial e and I profiles or to the boundary layer edge 
boundary conditions on e and 5 - as long as reasonable values 
are used. 

2.3.2 Boundary Conditions 

Since the inviscid forms of the four models are considered, it 
is not possible, to integrate through the viscous sublayer. 
Hence, it is necessary to assume that the sublayer has zero 
thickness and then match to the law of the wall. A singular 
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perturbation solution of the jQodel equations, valid as y^O, 
yields the proper boundary conditions. Such a perturbation 
solution (valid for boundary layers with and without pressure 
gradient) for all four models is listed below. 


u , 
e 


1 


== ..2 


u^/a* 

T 


Z = /a*' 


Ky 


( 69 ) 


The quantity u^ is the friction velocity, ic is Karman’s con- 
stant, a* is the constant in the Saffman formulation and B is 
the constant in the law of the wall. In all calculations, the 
constant B was held fixed at 5*5. 

The boundary conditions in Equation 69 must be applied close to 
the surface, viz, at least closer than a y"^ (defined as u^y/v) 
of 20 and preferably at a y"*" of about 10. Applying the condi- 
tions closer than y"*" of 10 made little difference while apply- 
ing them above 20 significantly affected the solutions. On 
first inspection, it may seem inaccurate to apply the boundary 
conditions so close to the surface that the law of the wall does 
not actually hold (i.e., y^£30). However, Equations 69 are 
thgorigorDususolutlontto theilnMscidembmentumaequationyas, y 
approaches 0. That is, by assuming that the sublayer has neg- 
ligible thickness. Equations 69 are the singular perturbation 
solutions for the various models and are the appropriate 
boundary conditions in the limit as y approaches 0. 


The perturbation analysis of the defect layer described in Sub- 
sections 2.1 and 2.2 indicates the appropriate boundary-layer- 
edge value of the turbulent mixing energy. The analysis shows 
that a value of e corresponding to a*e/u^ = .01 is the largest 
value consistent with the equations that yields accurate defect 


30 



layer velocity profiles; using smaller values of emake little 
difference. Since the larger value is somewhat eastier to handle 
numerically, all computations were performed with a*e/u^ = .01. 

Because the model equations predict existence of a sharp turbu- 
lent-nonturbulent interface the bpundary-layer edge, imposing 
95-/9y = P (Equation 35) proves to be very difficult. Hence, con- 
sistent with mixing-length approximations, a value of S,/j/o^6 = .09 
at the boundary layer edge was imposed. Varying the edge value 
of over the range of .06, to .12 made little difference 

in the numerical solutions. Hence, an important conclusion is 
that the models are not inordinately sensitive to the edge values 
of e and i. 

2.3*3 Initial Profiles 

Prior to this study, determining reliable initial profiles for 
the turbulence quantitldsbwas always an area of uncertainty. 

As will now be- demonstrated, this uncertainty has been eliminated. 

Starting profiles for all models are most conveniently formulated 
in terms of the length scale, 1; Equatlpns 6^ and 68 are then 
used to determine w and e^. Noting that e is interpreted as 
9/4<v’^>, starting profiles can be easily obtained if <v*^>, 
<-u’v’>, and u data are available at the starting location. 

That is, e can be determined directly from the <v’^> data while 
i can be found from the following equation 

5, = <-u’v’> (70) 

dy 

Howdver, such complete information is not always available. 
Nevertheless, good profiles can still be determined with only 
<!?iiifvJ >aanduudd&.feaaattlbhessfeartiinglloeati.on . TThisiis 'done using 
the following procedure. 
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First, it is assumed that an arbitrary e starting profile can 
be written as 




( 71 ) 


where ’i' is a constant to be determined and ©ppgL high- 

Reynolds-number PPBL e profile.^ Prom empirical observation, 
this is a good approximation for equilibrium boundary layers. 
The PPBL 9/4<v*^> profile can be approximated by 


'PPBL 


u 


a“ 


COS' 


(ID 


(72) 


Next it is noted that as y->-0, e^<u^y. Since e also equals 

any augmentation in the e profile must be reflected in 
the I profile as well. Therefore, e can be rewritten as 


e 


e 


1/2 


PPBL 



^PPBL 


so that we must approximate the H profile as 


Z 



Spbl 


(73) 


where high Reynolds number PPBL length scale which 

empirically has been found to have the following character: 


^PPBL 


y y £ • 096 /k 

.O9v/o^6 y > .096/k 


(74) 


Therefore it is only necessary to determine the value of 4' for 
each particular flow. It is known that across a PPBL, <-u’v'>/ 
9/4<v' ^ > w 0 . 3 • This implies that the <-u*v’> and <v.’^> profiles 
are similar. It is possible that this is nearly the case for 


^ Note that Equation 71 fails to satisfy the boundary condition 
e^u^/a* as y-^0.. Experience has shown that the computations correct 
this error in a reasonably short distance. 

J 
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other flows. If this is indeed true, then it is a straight- 
forward operation to find 'i'. If the <-u’v’> and <v’^> profiles 
are similar, then the ratio of the maximum value to the wall 
value is the same for each profile. This ratio has the value 
(see Figure 8). 

This procedure was studied in great detail for the three cases 
for which extensive data are available, viz, Andersen’s adverse 
pressure gradient flow, Bradshaw’s adverse pressure gradient 
flow, and PPBL flow. For example, in the Andersen case, only 
<-u’v’> data are available at the initial location, x=2 ft' (.61 m) 
However, <-u ' v ' >//9/^<v’ « 0. 26 almost uniformly for y/6 < 3/^ at 
x=6 ft ( 1.83 m) . Hence, the profiles were assumed nearly similar 
at the starting location, x=2 ft (.61 m). The ratio of the maxi- 
mum value of <-u’v’> to the surface shear stress at x=2 ft (.61 m) 
is 1.5- The e and i£ profiles were then generated using = 1.5* 

As a check on this procedure, we noted that Andersen also lists 
mixing-length profiles at x=2 ft (.61 m) . The Inferred values of 
i are compared with the mixing-length data in Figure 9; agreement 
is excellent. 

Sufficient data are available to directly obtain e and £ profiles 
for the Bradshaw adverse pressure gradient flow. Using 4'=1 gives 
a length scale profile close to the Bradshaw data with the excep- 
tion of the measured decrease in £ near the boundary layer edge 
as shown in Figure 10. The Bradshaw case was run both with the 
actual £ profile and with the profile constructed with 'l'=l; 
results differed only slightly. Figure 11 compares the Initial e 
profile constructed with f=l with experimental data; the figure 
shows that the cos^ fit is quite accurate. 

The £ profile used in the FPBL computations is shown in Figure 12. 
This profile does not quite coincide with the one constructed with 
'}'=1. The difference occurs because the FPBL flow was initiated 
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Figure 8. Use of similar <-u'v’> and e profiles 
to determine H'. 
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y/6 


Figure 9. Comparison of Inferred and measured Initial 

mixing length profiles for the Andersen adverse 
pressure gradient flow: x=2 ft (.61 m). 





Figure 10 . 


Initial length scale profile for the Bradshaw 
adverse-pressure-gradient flow; x=2 fteC .61 in 


). 
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Figure 11. Initial turbulent mixing energy profile 

for the Bradshaw adverse-pressure-gradient 
flow: x=2 ft-eC.6l m) . 






at a very low Reynolds number; hence the peak length scale had 
to be enhanced to be consistent with empirical observations for 
low Reynolds number flows. The FPBL cases were also computed 
with a profile using 'l'=l; a transient was present in the early 
part of the computation but quickly settled out. The standard 
' EPBL e profile defined in Equation 72 was used in both runs. 

In summary, by using both the inferred and actual Z profiles, 
we have shown that the computations are not overly sensitive 
to the starting profiles as long as a reasonable peak length 
scale is used. Hence, the starting profiles defined by Equa- 
tions 71-74 are suitable for equilibrium boundary layers. 

Simple PPBL profiles ('l'=l) were used for both the Ludwleg- 
Tlllmann and So-Mellor flows. 

2.3*4 Computations 

As noted earlierpi-ncthls.gSubseGbiQnvhthgaSafjfiman^ic-BaSeline , Ng- ■ 
Spalding, and Jones-Launder turbulence models were applied to 
the following five flows: (1) FPBL flow; (2) Bradshaw adverse 

pressure gradient flow; (3) Andersen adverse pressure gradient 
flow; (4) Ludwleg-Tillmann favorable pressure gradient flow; 

(5) So-Mellor constant pressure flow over a convex wall. The 
results of these computations are shown in Figures 13-20. 

All of the models are accurate for FPBL flow. Computed proper- 
ties for the Baseline, Mg-Spaldlng, and Jones-Launder models 
agree closely with measurements, while Saffman-model predictions 
generally differ from corresponding data by about 10%. The 
Baseline model yields a velocity profile in linear coordinates 
which agrees most closely with experimental data (see Figure 13)* 

The Baseline and Ng-Spalding models yield the closest agreement 
between computed and measured properties for both adverse pres- 
sure gradient flows, whereas Jones-Launder-model and Saffman-model 
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lUOUYE CORHELATION 


WIEGHARDT DATA 


Skin Prlction 


Shape Factor . - ^■• 


KLEBANOFK DATA 
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X 
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' . <'v:/r^r > i- ; 

Velpcity oE^pf lie ;; 


Figure 13 . 


Comparison of computed and measured flow properties 

for flat-plate boundary-layer flow; '’ Baseline, 

— ‘ — Saffman, Jones-Launder and Ng-:Spaldlng. 





COLES ’ VERSION OP BRADSHAW 
DATA 

BRADSHAW DATA * 


BRADSHAW DATA 


Skin Friction 




Shape Factor 


- •<> ;> ' .'i- Iv ’ ■ 




BRADSHAW DATA 


/ I I 

■;; ' 1-35 




^ EHADSHAW DATA 

A COLES’ VERSION 0? BRADSHAW DATA 


§ /T' 
V H 




Velocity Profile j x=7 feet 


Velocity Prof lie; x=7 Peet 


Figure 14. 


Comparison of computed and measured flow properties 
for the Bradshaw adverse-pressure-gradlent flow; 

•VJiaseli-he§,e4-T • -trr, Saffman.yn , .-, Ng-Spaldlng, 

Jones-Launder .' 







PfgjireelS. (af!pcbJ^mfee.dneanidraBfe£sjar.ed flow properties 

for the Andersen adverse-pressure-gradient flowj; 

r-r^ WBaseilSe^clr^ Saffman, — -- Ng-Spalding, 

Jones-Laiinder . 
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> Velocity Profile ; 3^14 .4 feet-, 


Figure l6. Comparison of computed and measured flow properties 
for the Ludwleg-Tillmann favorable-pressure- 

gradient flowj B'aseldne'.',.^r-r-.-^ — Saffman, 

— Ng-Spaldlng, Jones-Launder . 
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Piggre-il'T'- C^cCipmp’a-ciS'STicofpvdoinputed and measured (e = 9/4<v'^>) 
MrPuieHfce®txing^energx -profiles ; — . Baseline, 

^ — Saffman, Ng-Spaldlng, ^ — Jones- 

Launder. 
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Fi[gSi£?:eel'8 . Comparison of computed and measured length 

scale profiles;; ^iBh-'S^ellfte,' 

Saffman, Ng-Spaldlng, 

" Jones-Launder . 
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predictions differ from the data substantially. Most particu- 
larly, note the excellent agreement between experimental and 
computed velocity profiles in linear coordinates for the 
Baseline model (see Figures l4 and 15)* The Ng-Spaldlng-model 
velocity profiles in linear coordinates do not agree nearly as 
well with the data. However, with the exception of the veloc- 
ity profile in linear coordinates, the Ng-Spalding model yields 
closest overall agreement to the Andersen data. 

For the Ludwieg-Tlllmann favorable pressure gradient flow, the 
Baseline, Jones-Launder and Ng-Spalding models yield predic- 
tions very close to the experimental data (Figure 16). The 
Saffman model does not. do quite as well; the velocity profile 
in linear coordinates serves as a reminder of the Saffman model's 
deficiency in the defect layer. 

Figures 17 and l8 compare computed e and i profiles with experl-' 
mental data for 9/4<v.*^> and mixing length, respectively. 
Consistent with the defect-layer analysis of Subsections 2.1 
and 2.2, in all three cases shown the Baseline model has length 
scale profiles much closer to the measured mixing-length pro- 
files than does the Saffman model. All models do a good Job of 
reproducing the 9/4<v'^> data indicating that the efficacy of 
a two-equation model is controlled by the accuracy with which 
the turbulent length scale can be computed. 

Of these four flows, the models do poorest on the Andersen flow. 
However, the Andersen flow has a relatively low Reynolds number 
and, without suitable viscous corrections, the models are realis- 
tically only applicable at high Reynolds numbers. This is best 
exemplified in Figure 15* The Baseline and Ng-Spalding models 


i* 1 

The Saf fman-model solutions generated by Wilcox and Chambers 

(see Figure 1) were done with a value for Zq/Z^S of less than 

.02; although accurate skin friction and shape factor were 

obtained, the inflected profiles shown in Figure 1 attended the 

very small value of Jig. 
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stray from skin friction data early in the computation. How- 
ever, agreement is much better by the end of the computation, 
where the boundary layer is finally approaching a well-developed 
turbulent state. 

As a final part of the comparisons, the models* ability to pre- 
diet effects of streamline curvature was tested. Computed skin 
friction for So-Mellor's constant pressure flow over a convex 
wall (Figure 19) shows that two-equation models in their orig- 
inal form fail for flows with streamline curvature. However, 
as shown by. Wilcox and Chambers,^ better agreement between theory 
and experiment can be obtained by adding a curvature correction 
term to the turbulent* mix j ng-energy equation. Using this curva- 
ture term, the So-Mellor flow was recomputed with the Baseline 
model. As shown in Figures 19 and 20, close agreement between 
calculated and measured skin friction and velocity profiles is 

f ' 

then obtained. Thd&cbrreetionnterm'nwill ■be'discus'-Sed''^ln- detail 
in Section 3 when other flows with streamline curvature are 
considered. 

2.3*5 Summary 

Fbrrthe flows without streamline curvature, the Baseline and 
the Ng-Spaldlng models yield results in much closer agreement 
with the data than do the Jones-Launder and Saffman models. 
Furthermore, the Baseline model's velocity profiles in linear 
coordinates are in very close accord with the data and consis- 
tently closer than any of the other models'. Hence, the 
Baseline model appears to be the most accurate of the various 
two-equation models. Computation of the curved-wall case shows 
that none of the models in their original form accurately pre- 
dicts effects of streamline curvature. However, by adding a 
streamline curvature correction term, accurate predictions can 
be obtained with the Baseline model for flow over a convex 
wall. This curvature correction will be shown to work well for 
a number of flows with curved streamlines in Section 3- 
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2.4 


VISCOUS COMPUTATIONS 


The utility of the Baseline model, would be severely limited if 
it were applicable only where the law of the wall holds close 
to a boundary. It would be inapplicable for transitional and 
separating flows for example and, more generally, it would per- 
form poorly for boundary layers with arclength Reynolds numbers 
less than a million. Hence, encouraged by the success of the 
Baseline model in the matching computations, viscous calcula- 
tions were performed. The objective of performing viscous com- 
putations is to test the model's ability to predict low Reynolds 
number effects. The four flows with no streamline curvature con- 
sidered in the preceding Subsection were recomputed. In addition, 
a viscous solution of a flat plate boundary layer with uniform 
blowing was obtained. Table 2 lists the specific flows con- 
sidered and the corresponding data sources. 


Table 2. Incompressible Viscous Calculations 




.PLOW 


DATA:.SaUR,CE,S:s 


Plat plate boundary layer 


• Bradshaw adverse pressure 
gradient 

• Andersen adverse pressure 
gradient 


Same as in matching 
calculations ■ (Table 1 ) . 


• Ludwieg-Tlllmann favorable 
pressure gradient 

• Andersen flat plate with 
uniform blowing 


Andersen 


4 


data. 


In order to proceed with viscous computations, viscous modifica- 
tions to the model equations had to be Included. Wilcox and 
17 

Traci used perturbation techniques to study the sublayer of 
a turbulent boundary layer. They found that in addition to adding 
molecular dlffuslbp to the momentum and model equations, the fol- 
lowing straightforward modifications to a and a* are needed to 
obtain accurate viscous computations: 
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( 75 ) 


2|= l.(l.X)e^%4 

oo 

JL = 1- (I-X) e-”®T/2 
a_ 


where 

X = 1/11 , a*= 3/10, a„= 3 /a* - 2a<^ = .'3319 


and Re^p is the turbulent Reynolds number defined by 



(76) 


Boundary-layer-edge boundary conditions and initial profiles were 
the same as those used in the preceding subsection; for the blow- 
ing case, FPBL initial profiles were used. At the surface, the 
no-sllp boundary condition was. imposed while e and to were given 
by 



(77) 


The quantity S is a function of surface roughness and blowing 

-| Q 

lyelocity and, as shown by Traci, 


where 


s-^ 




(78) 


^B 


1 + v^/u^ 


(79) 


51 



and 



' 36 

ku /v 

^ T 


, 1/2 


ku/v 


C80) 


In Equations 79 and 80, is blowing velocity and k is roughness 
height. The nonblowing computations were performed with very 
small values of k to simulate a perfectly smooth wall; we have 
found empirically that the smooth wall limit in which is 

very closely approximated with finite values of S in excess of 
about 300 - 

Figures 21 through 27 show the results of the viscous computa- 
tions for the boundary layers with no mass injection. The vis- 
cous solutions for the FPBL, Bradshaw adverse pressure gradient, 
Andersen adverse pressure gradient, and Ludwleg-Tillmann flows 
are very similar to the Inviscld Solutions (Figures 13-18). For 
the blowing case, skin friction corresponds well with data, com- 
puted shape factor is within 10 % of data, and velocity profiles 
are within 15 % of the data (Figure 25). 


Results of the FPBL and Andersen adverse pressure gradient flows 
best exhibit the model's ability to accurately predict viscous 
effects. For FPBL flow. Figure 21 shows close agreement between 
the predicted and experimental shape factor for a wide range of 
Reynolds numbers. Additionally, the figure exhibits the Baseline 
model's ability to handle viscous effects by displaying velocity 
profiles in sublayer coordinates for three different Reynolds 
numbers. Note that the Andersen adverse pressure gradient calcu- 
lation was initiated near the leading edge of the plate, whereas 
the matching calculation was begun at a plate length of 2 ft (.61 m) 
By including viscosity and integrating from near the leading edge, 
the results are in much closer agreement with measurements than 
those obtained by matching to the law of the wall. Figure 23 
shows velocity profiles in sublayer coordinates for two different 
Reynolds numbers. 
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Velocity Profile; x=7 feet 


Velocity Profile; x=7 feet 


Figure 22... Comparison of computed and measured flow properties 
for the Bradshaw adverse pressure gradient flow. 
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Velocity Profile; x=6 feet 




. Velocity i^of lies 


Figure 23. Comparison of computed and measured flow properties 
for the Andersen adverse pressure gradient flow. 
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.5 


LUDWIEB AND TIIXHANN DATA 

COLES' VERSION OP LUDWIEO AND 
TILLKAKN DATA 


LUDWIEO AND TILLMANN DATA 


6 8 


Skin Friction 


Shape Factor 




0 .2 .6 .8 1.0 


Velocity; :pi’ofile .x=.l4.4. ‘feet ee; 


Velocity jRrof ile ; x=l4 . 4 feet 


Figure 24. Comparison of computed and measured flow properties 
for the Ludwieg-Tillmann favorable pressure grad- 
ient flow. 






Velocity profile; x=f0, Inches Velocity Profile; x=90 Inches 


Figure 25 . Comparison of computed and measured flow properties 
for flat-plate boundary-layer flow with uniform 
blowing. 
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COMPUTED 


COMPUTED 

• BRADSHAW DATA 



y/6 



, Andersen Adverse^ Vp Flowj. jc=6 feet 


}B%g,ure?626 .CcffiSo,rat);a'rilsf 6 n''Of cdmputedsand- measured (e = 9 / 4 <v’^>) 
turbulent mixing-energy profiles. 





& 

y/a^6 


I 




PPBL Plow; Re =^.20*10® Bradshaw Adverse- Vp Plow; x=7 feet 

X 

£ 

VoS’s 



Andersen Adverse-Vp Flow; x=6 feet 


Figure 27. Comparison of computed and measured 
length scale profiles for flat-plate 
boundary- layer flow. 
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Results, using the. Baseline .model, for both the inviscid and 
viscous calculations of equllibr.iLun boundary layers are most 
encouraging. The inviscid results, indicate, that the Baseline 
model performs at least as well as any of the other models and 
usually better. The viscous computations show that the model 
equations can easily and successfully be. integrated through 
the viscous sublayer and that they yield solutions that corre- 
spond well with data for flows with adverse, zero, and favor- 
able pressure gradients and for flows with blowing. Since the 
Baseline model has these capabilities, it can be regarded as a 
reliable foundation for analyzing flows with streamline 
curvature and coordinate system rotation. 


60 



3. BOUNDARY LAYERS ON CURVED SURFACES 

We now turn our attention to effects of streamline curvature 
and system notation on turbulent boundary layers. As stated 
in Subsection 2.3 - the Baseline model in its original form 
does not accurately predict the effects, of streamline curva- 
ture. However with the addition of a straightforward curva- 
ture correction term, the Baseline model can successfully be 
applied to flows with curved streamlines'. This section gives 
a brief overview of the development of the correction term and 
tests It for boundary layers with curved streamlines and for 
rotating channel flow. Both compressible and Incompressible 
flows are considered. 

3.1 STREAMLINE CURVATURE/SYSTEM-ROTATION MODIFICATIONS 

With an understanding of some of the' physics of flow over 
curved walls and the proper Interpretation of the turbulence 
quantity, e, a straightforward curvature correction can be 
derived. This is explained In detail by Wilcox and Chambers.^ 
Hence only a brief description will be presented here. 

Three points are essential to the argument. The first impor- 
tant point is that the mixing energy is directly proportional 
to <v’^> (the component of the fluctuating velocity normal to 
the surface); m is the rate at which <v’^> kinetic energy is 
converted to other modes of energy (e.g., mean kinetic energy, 
thermal energy, other fluctuation modes). The second key point 
is that the equation for e should be regarded as the <v’^> 
component of the full Reynolds stress equation^. Tine- -r:f Inal — 

key observation is that if the tensor transformation properties 
of the <v'^> equation are endowed upon the equation for e, 
straightforward curvature (and system rotation) correction 
terms can be added which account for altered turbulence structure 
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Implementing Point 3 proceeds as follows. In a curved-stream- 
line coordinate system with local radius of curvature, R, the 
equation for the Instantaneous vertical velocity, v, has a 
Coriolis acceleration term, viz. 


If + u-W - f = 


+ diffusion 

P 9y 


( 81 ) 


Hence, multiplying Equation ,83'^by v' and time averaging, we 
obtain 


+ u*V<v.’^/2> - ^<u’v'> = other terms (82) 

at — n 

where, for simplicity, u, denotes mean value. Then, defining 
e E 9/4<v' ^>. and approximating <-u’v’> = — Equation 82 
becomes (Cw^lt'h’&,overhea^d*-.:bars>-d.l?oppeidpf-c^^ :6o'Hven4'ehee'!i)> c*' c J 

If + u-Ve = -flllf- + other terms (83) 

The term multiplied by 9/2 is the curvature correction term. 
Similarly, for flow in a coordinate frame rotating with angular 
velocity we obtain 

+ u*Ve ==— 9^^— -1^ + other terms (84) 

ot — ' CO dy 

The e and co equations used for the boundary layer flows with 
curved streamlines and for rotating channel flow are hence 
assumed to be: 


Note that the correction term was incorrectly written by 
Wilcox and Chambers^ as being proportional to 9- 



Boundary Layers with Curved Streamlines 


Turbulent M ixing Energ y 


9© I 9e p « I 9u 

^93E 


9y 


9y 


u 

R 


- «*“> ^ - I I S It 


Dissipation Rate 


u 


^ -If = t“l|t - gl - ^ ^ 


Rotating Channel Flow 


Turb ule n t Mixing En e r gy 


dy 0 ) dy 


(a*lf^l - 3*w)e 


9 ^ = 

(jo dy 


DlsslpaW Hatj 

f[(vfa|) + [a|M| . , 3 


3.2 INCOMPRESSIBLE APPLICATIONS 

The Baseline model with the curvature correction term was first 
applied to several Incompressible flows (Table J). 


Table 3- Incompressible Viscous Calculations 
for Boundary Layers with Curvature 


PLOW 

DATA SOURCE 

• Low Reynolds number turbulent 

19 

Patel data. 

flow past a cylinder 


• Constant pressure flow over a 

So-Mellor^^ data; Meroney^^ 

convex wall 

correlation. 

• Separating flow over a E<. 

)-S 0 ^Mellor^^ data. 

convex wall 


• Rotating channel flow 

Johnston^*^ data. 


(85) 


(. 86 ) 


(87) 


( 88 ) 
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The most difficult test of all the Incompressible cases 
considered is low Reynolds number turbulent flow past a cylinder. 
In addition to having streamline curvature and adverse pressure 
gradient, the flow is dominated by viscous phenomena. Results 
of the computation are shown in Figures 28-30. All computed 

quantltltes are in close agreement with the experimental data . 

19 

of Patel. Computed and measured velocity profiles are com- 
pared in Figure 28. The largest discrepancies are less than 
10^ throughout the flow. Even the predicted and measured sep- 
aration points are reasonably close; the computed separation 
occurs at circumferential angle (t)=101° while separation was 
observed to occur at (j)=110°. As shown in Figures 29 and 30,. 
computed shape factor and momentum-thickness Reynolds number 
differ from corresponding measured values by no more than 6^, 
even as the flow nears separation. 

Figures 31-33 exhibit computed and measured flow properties for 
flow over a convex wall with constant pressure and with adverse 
pressure gradient. Figure 31 shows that computed and measured 
skin friction are within a few percent. In the separating flow, 
separation was observed between 5-50 and 5-83 feet (1.68 and 
1.78 m) downstream of the leading edge, while the computation 
indicates separation will occur a little farther downstream. As 
seen from Figure 32, predicted shape factors show somewhat larger 
discrepancies from the measured values; the largest discrepancy 
is 12^. Finally, Figure 33 shows that computed and measured 
velocity profiles differ by no more than 6/5. 

The final incompressible application is to rotating channel 
flow (Figures J|i-36). Figures 3^ and 35 show the contrast 
between nonrotating and rotating channel velocity profiles. 

Figure 36 presents the ratio of the rotating friction velocity, 
u^, to the nonrotating friction velocity u.^^, for a range of 
rotation numbers. The predictions are generally within 10% of 
corresponding measurements. 



6 IN DIAMETER 
CYLINDER 


8 CM 8 


OBSERVED 

. SEPARATION POINT 

= 101 " 1 

COMPUTED ^^^a-a-a^^_a-a.a^-a-a 
SEPARATION K 


0° 20° 40° 60° 80° 100° 120° 140° 160° l80° 

(j) 

Exper?3jneritaiL:ta35rangement :"Of- Gylinder-and ;measoi’edi;pEfessu3?eJ' distribution 

0 0.2 0.4 0.6 0.8 1.0 


60 ° 70 


A A PATEL DATA 
COMPUTED ' 


iwi 


Velocity profiles at several stations on the cylinder. 


Figure :2.8. 


Comparison of computed and measured flow properties 
for Patel's low-Reynolds-humber turbulent flow 
past a cylinder; Re^ = 5 • 01 • 10 ® . 
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H 



Figure 29. Comparison of computed and measured 

shape factor for Patel’s low-Reynolds 
number turbulent flow past a cylinder 

= 5 . 01 *10 
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X (ft) 


Figure 31* Comparison of computed and measured skin 
friction for flow over a convex surface. 
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H 



X (ft) 


Figure S?. Comparison of computed and measured 
shape factor for flow over a convex 
surface . 
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u/u 
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Figure Comparison of computed and measured, 

velocity profiles for flow over a 
convex surface. 
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u/U 


m 



y/H 


Figure 34. Comparison between computed and measured 
flow properties for non-rotating channel 
flow. where Uj^ is the average 

flow velocity. 
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Figure I5. Comparison between computed and measured 
flow properties for rotating channel 
flow. R = where U Is the average 

flow velocity. 





JOHNSTON 

■ .23 -,¥0iD 

• . Re = 25,700 
A Re = 33,100 
4 Re = 11,400 

COMPUTED 


Figure 36 . Comparison of computed and experimentally 
measured values of u^-Zu^- as a function 
of the rotation number, ° fiH/Ujjj f or rotat- 
ing channel flow; = rotation speed, 

H = channel width, Ujjj = average velocity. 

Re = Reynolds number based on the average 
velocity . 
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3.3 COMPRESSIBLE APPLICATIONS 

As a final test of the theory, the Baseline model was applied 
to two compressible flows (Table 4). Following Wilcox and 

Table 4. Compressible Viscous Calculations 
with Streamline Curvature 


FLOW 

DATA SOURCE 

• Mac ]3 2, ,5 flow witfe heat 
transfer over a convex wall 

• Mach 3.5 nearly adiabatic 
flow over a concave wall 

22 

Thomann data. 

23 

Sturek-Danberg data. 


21 

Alber, we Introduce the specific turbulent dissipation rate, 
fi, so that the boundary layer form of the compressible model 
equations with the curvature modification are as follows. 

Specific Turbulent Mixing Energy 


pu 


3e 

3x 


+ pv,|f.= [o.*||f-^|-8*pa]pe+^C(. 


+ 


96 -1 ^ ^ ^ 9u 
ay-* 2QR 9y 


Specific Turbul.erit Dissipation Rate 


]pS!)p!2- 


9 r ( X ^ ^ ^ ^ “I 

ay^^^ ^ 3y ^ 


The quantities p and p denote density and molecular viscosity 
(given by Sutherland’s law), respectively. 

The Baseline model was first applied to the Thomann flow using 
a wall boundary condition on $1 that simulates a smooth wall 
(the case of actual Interest). As a generalization of the 


(89) 


(90). 
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Incompressible boundary condition (Equation 77) j the wall 
value of the dissipation rate is given by 


Q = 

w 


Su2 


( 91 ) 


where S is the universal function defined in Equations 78-80 . 
and is the molecular viscosity at the wall. Using a value 
of S corresponding to a smooth wall (S of the order of 10**) 
gave skin friction that was below the Hopkins-Inouye compres- 
sible skin friction correlation by about 20^. Also, for this 
value of S, the Stanton number, S^, was as much as lS% too low 
^>Figure 37). However, it was found that by simulating a rough 
wall (S on the order of 10), good skin friction and Stanton 
number distributions for the equilibrium portion of the boundary 
layer were obtained as shown in Figures 37 and 38. 


The necessity of simulating flow over a rough wall is not as 
yet understood and is an area for future research. However, 
by at least assuring good equilibrium solutions, it was pos- 
sible to examine how the model behaves when the boundary layer 
is disturbed by curvature and/or pressure gradient. 

In the region where curvature begins, consistent with measure- 
ments, the model predicts a 10% reduction in Stanton number 
(relative to the plane-wall value). However, the model fails 
to predict the observed rapid decrease in beyond x=1.25 ft 
(.38 m) . The discrepancy may be due to the fact that curva- 
ture affects the heat transfer differently from the way it 
affects skin friction. 


The Sturek-Danberg calculation proceeded as did the Thomann 
case. The appropriate large value of S resulted in skin fric- 
tion as much as 19% iliower than the Hopkins-Inouye skin friction 
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f 



Re 

X 


Figure 38 . Comparison of computed skin friction with 

the Hopklns-Inouye correlation for Mach 2.5 
flow over a convex surface. 
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correlation (Figure 39 ). Hence, as in the Thomann case, the 
computation was performed with S=10. As shown in Figure 39 , 
the predicted c^ is very close to the Hopkins-Inouye value, yet 
is considerably higher than the measured values in the equil- 
ibrium portion of the flow. However, the data are suspect in 
this region since even Sturek and Danberg note that their 
measuring techniques were only accurate to within 15^. 

The predicted overall decrease in surface shear stress is 
consistent with the measurements although differences of about 
20%, of scale are present. The discrepancies in the curved 
region of the flow may also be caused by the inaccuracies in 
the measurements so that this case does not provide a definitive 
test of the theory. 
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Figure 39. Comparison of computed and measured skin friction for 
Mach' 3.5 flow over a concave surface. 




4. DISCUSSION 


Perhaps the most important conclusion that can he drawn is that 
the Baseline model is at a state of dg^velopment where it can 
be used for general engineering applications. In particular, 
results of the viscous equillbrlum-boundary-layer computations 
in Section 2 show that the Baseline model is as accurate as 
mixing-length theory where iriixlng- length theory is expected to 
apply- Furthermore, because of the model's generality and 
more-solid physical foundation, its range of applicability 
extends far beyond that of mixing-length theory's. The model's 
range of applicability includes flows with curved streamlines 
and coordinate system rotation. Using the Baseline model, 
excellent quantitative agreement between theory and experiment 
has been obtained for effects of curvature and rotation on 
incompressible boundary layers. Although more work remains to 
be done for compressible boundary layer applications, results 
for the two supersonic flows considered in Section 3 are also 
encouraging. 

Consistent with the project objective, the path is now cleared 
for proceeding to analysis of streamline-curvature"''''and ' system- 
rotation effects on turbulent boundary layers. Presumably, 

RSE model deficiencies can now be eliminated by changes analo- 
gous to those made to the Saffman model which resulted in the 
Baseline model. Such changes should be made. Work with the 
Baseline model should continue as well. That is, parallel 
development of the RSE and Baseline models appears to hold the 
most promise for yielding useful, accurate engineering tools 
for predicting curvature/rotation effects. 
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APPENDIX 


As noted In the Introduction, numerical solutions with the RSE 

model were generated in this study. The model was applied to 

'24 

Bissonnette and Melloris ffl'OWppast -a segmented .cylinder. "In- 
the Bissonnette-Mellor experiments, the forebody is stationary 
and the aft body rotates at a constant angular velocity. The 


fact that the pressure is not constant in the direction normal 

to the surface (- 1 ^ 7 ^ 0) was taken into account in the computations 

dy 



The computed mean flow properties and turbulence quantities are 
shown in Figures A1-A7; experimental data of Bissonnette and 
Mellor are included for comparison. As shown in Figure Al, pre- 
dicted streamwlse (6f ) and azimuthal (cf ) skin-friction 

j- X z 

coefficients are both about 20^ higher than measured. Similarly, 
the various computed displacement and momentum thicknesses are 
generally 20^ higher than measured (Figures A2 and A3). Fig- 
ures A4-A7 present velocity and Reynolds stresses on the swirling 
afterbody. The figures show poor agreement between experimental 
and computed properties in almost every case. Inspection of the 
streamwlse velocity profile (Figure A4) shows that the worst 
agreement is in the defect layer of the boundary layer. This 
is consistent with the results obtained in analyzing Saffman- 
model-predlcted defect layer structure (see Subsection 2.1). 

That is, examination of Figure a 4 shows that the computed 
streamwlse velocity gradient, 9u/9y, is much larger than mea- 
sured in the defect layer. Too strong a gradient results in 
too much production in the various Reynolds stress equations, 
and as a result the Reynolds stresses are overestimated. With 
the exception of <-u’w’>, the Reynolds stress components are 
overpredicted as can be seen in Figures A6 and A7. 
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The large discrepancies between theory and experiment for the 
rate of growth of the swirling boundary layer are directly 
attributable to the excess production of turbulent energy. 

High turbulent energy levels are synonymous with greater tur- 
bulent mixing; because of the overpredicted mixing, the swirl- 
ing boundary layer grows too quickly and hence the swirling 
displacement thickness, 5*, and swirling angular momentum 

thickness, 0 , are too large, 

xz 

In summary, results of this computation lend further credence 
to the claim that the problems plaguing the RSE formulation 
have the same origin as the deficiencies of the Saffman model. 
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Swirl turned on 


Computed 
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= ''x/1/2 pu* 

= T /1/2 pu^ 

^w ® 
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Figure A4. Comparison between computed and measured 
streamwlse velocity profile at x=3.17 ft 
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